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and where it decays exclusively into a charm quark and a neutralino, 95% confidence-level upper 
limits are obtained in the scalar-charm-neutralino mass plane such that, for neutralino masses below 
200 GeV, scalar-charm masses up to 490 GeV are excluded. 
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The results of a dedicated search for pair production of scalar partners of charm quarks are 
reported. The search is based on an integrated luminosity of 20.3 fb~^ of pp collisions at ^ = 8 TeV 
recorded with the ATLAS detector at the LHC. The search is performed using events with large 
missing transverse momentum and at least two jets, where the two leading jets are each tagged as 
originating from c quarks. Events containing isolated electrons or muons are vetoed. In an A-parity- 
conserving minimal supersymmetric scenario in which a single scalar-charm state is kinematically 
accessible, and where it decays exclusively into a charm quark and a neutralino, 95% confidence- 
level upper limits are obtained in the scalar-charm-neutralino mass plane such that, for neutralino 
masses below 200 GeV, scalar-charm masses up to 490 GeV are excluded. 


PACS numbers: 14.80.Ly,13.85.Rm 

Supersymmetry (SUSY) [1-9] is a theory that extends 
the Standard Model (SM) and naturally resolves the hi¬ 
erarchy problem by introducing supersymmetric partners 
of the known bosons and fermions. In the framework of a 
generic i?-parity-conserving minimal supersymmetric ex¬ 
tension of the SM, the MSSM [10-14], SUSY particles 
are produced in pairs and the lightest supersymmetric 
particle (LSP) is stable, providing a possible candidate 
for dark matter. In a large variety of models, the LSP is 
the lightest neutralino, Xi. 

The scalar partners (squarks) of various flavors of 
quarks may, rather generally, have different masses de¬ 
spite constraints on quark flavor mixing [15]. Recent 
searches disfavor low-mass top squarks (stops), sbottoms, 
and gluinos, so direct scalar-charm (c) pair production 
could be the only squark production process accessible at 
the LHC. Searches for c states provide not only a possible 
supersymmetry discovery mode but also the potential to 
probe the flavor structure of the underlying theory. 

Since no dedicated search for c has previously been 
performed, the best existing lower limits on c masses are 
obtained from searches for generic squark and gluino pro¬ 
duction at the LHC [16, 17], and from the reinterpreta¬ 
tion of LHC searches [18] for direct pair production of 
the scalar partner of the top quark followed by decays 
ii —)■ c -I- Xi. The top squark searches have a final state 
similar to that expected for scalar charm quarks, but 
are optimized for small — m^o mass differences, and 
so have good sensitivity to the scalar charm quark only 
when rric — m^o < mw- 

In this Letter, a dedicated search for direct c pair pro¬ 
duction is presented. The scalar charm quark is assumed 
to decay dominantly or exclusively via c —>■ c -|- Xi. The 
expected signal is therefore characterized by the pres¬ 
ence of two jets originating from the hadronization of the 
c quarks, accompanied by missing transverse momentum 
^^miss) j-egTjiiting from the undetected neutralinos. 

The ATLAS detector is described in detail else¬ 
where [19]. This search uses pp collision data at a center- 


of-mass energy of 8 TeV recorded during 2012 at the LHC. 
After the application of beam, detector and data quality 
requirements, the data set corresponds to a total inte¬ 
grated luminosity of 20.3 fb“^ with a 2.8% uncertainty, 
using the methods of Ref. [20] . 

The data are selected with a three-level trigger sys¬ 
tem that required a high transverse momentum (px) jet 
and A™*®® [21]. While events containing charged lep¬ 
tons (electrons or muons) in the search region are ve¬ 
toed, single-lepton triggers are used for control regions. 
Events are required to have a reconstructed primary ver¬ 
tex consistent with the beam positions, and to meet basic 
quality criteria that suppress detector noise and noncolli¬ 
sion backgrounds [22]. Jets are reconstructed from three- 
dimensional topological calorimeter energy clusters by 
using the anti-A:* jet algorithm [23, 24] with a radius pa¬ 
rameter of 0.4. The measured jet energy is corrected for 
inhomogeneities and for the noncompensating response 
of the calorimeter by using px- and ? 7 -dependent [25] cor¬ 
rection factors [26]. The impact of multiple overlapping 
pp interactions (pileup) is accounted for using a tech¬ 
nique, based on jet areas, that provides an event-by-event 
and jet-by-jet correction [27]. Only jet candidates with 
px > 20 GeV within \r]\ < 2.8 are retained. 

Electron candidates are required to have px > 7 GeV, 
|ry| < 2.47 and to satisfy “medium” selection criteria [28]. 
Muon candidates are required to have px > 6 GeV, 
|? 7 | < 2.4 and are identified by matching an extrapo¬ 
lated inner-detector track to one or more track segments 
in the muon spectrometer [29]. When defining lepton 
control regions, muons and electrons must meet addi¬ 
tional “tight” selection criteria [29, 30], and must satisfy 
track and calorimeter isolation criteria similar to those 
in Ref. [31]. 

Following this object reconstruction, overlaps between 
jet candidates and electrons or muons are resolved. Any 
jet within a distance AR = 0.2 of a medium quality 
electron candidate is discarded. Any remaining lepton 
within AR = 0.4 of a jet is discarded. Remaining muons 
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must have longitudinal and transverse impact parame¬ 
ters within 1 mm and 0.2 mm of the primary vertex, 
respectively. 

The calculation of is based on the vector sum 

of the calibrated px of reconstructed jets (with px > 
20 GeV and Ip] < 4.5), electrons, muons and photons, 
and the calorimeter energy clusters not belonging to these 
reconstructed objects [32]. 

Jets containing c-flavored hadrons without 5-flavored 
parent hadrons are identified using an algorithm, op¬ 
timized for charm tagging, based on a neural network 
that exploits both impact parameter and secondary ver¬ 
tex information and with a, B to D decay chain vertex 
fitter [33]. This algorithm achieves a tagging efficiency of 
19% (13%, 0.5%) for c-jets (6-jets, light-flavor or gluon 
jets) in tt events. The efficiency for tagging 6-jets is de¬ 
termined from measurements of dileptonic ti events [34]. 
The c-jet tagging efficiency and its uncertainty have been 
calibrated in inclusive jet events over a range of px using 
jets from collision data containing D* mesons [35]. Jets 
can be c-tagged only within the acceptance of the inner 
detector (jpl < 2.5), so only such central jets are retained 
after the above selection. 

Events are then required to have if™*®® > 150 GeV and 
one jet with px > 130 GeV to ensure full trigger effi¬ 
ciency, as well as a second jet with px > 100 GeV. The 
two highest-px jets are required to be c tagged. The mul¬ 
tijet background contribution with large caused by 

mismeasurement of jet energies in the calorimeters or by 
neutrino production in heavy-quark decays, is suppressed 
by requiring a minimum azimuthal separation (A^min) of 
0.4 between the direction and any of the three lead¬ 
ing jets. To reduce the effect of pileup, the third jet is 
exempted from this requirement if it has px < 50 GeV, 
IpI < 2.4 and less than half of the sum of its track px 
is associated with tracks matched to the primary vertex. 
In addition, the ratio of to the scalar sum of the 

transverse momenta of the two leading jets is required 
to be above one-third. Events containing residual elec¬ 
tron or muon candidates are vetoed in order to reduce 
electroweak backgrounds. 

After these requirements, the main SM processes con¬ 
tributing to the background are top quark pair and single 
top production, together referred to as top production, as 
well as associated production of WjZ bosons with light- 
and heavy-flavor jets, referred to as IV-l-jets and Z-|-jets. 
A selection based on the boost-corrected contransverse 
mass m-cx [36] is employed to further discriminate scalar- 
charm pair from top production. For two identical decays 
of heavy particles into two visible particles vi and V 2 , and 
into invisible particles, the contransverse mass [37] is de¬ 
fined as {[Ex(ni) -b Ex(w 2 )]^ - [Pt(i'i) - ■ 

The boost correction preserves the expected endpoint in 
the distribution against boosts caused by initial-state ra¬ 
diation. In the case of scalar-charm pair production with 
c —>■ c -b X]*, mcx is expected to have an endpoint at 


(to| — m?o)/we. For tt production, if both 6-jets are 
mistagged as c-jets, the mcx built using those two jets 
is expected to have a kinematic endpoint at 135 GeV. 

To maximize the sensitivity across the c-x5 mass 
plane, three overlapping signal regions (SR) are de¬ 
fined: mcx > 150, 200, and 250 GeV. The remaining 
tt background after the mcx requirement mostly com¬ 
prises events with one true c-jet from a W decay and a 
mistagged 6-jet from a top quark decay. Events in which 
a Z boson is produced in association with heavy-flavor 
jets where the Z boson decays into I'D also enter the high- 
mcx regions. The heavy-flavor jets often originate from 
a gluon splitting, g —>■ cc, which can lead to a small angu¬ 
lar separation between the resulting c-jets and therefore a 
small invariant mass rricc- The remaining tt background 
is also concentrated at low rricc- Gonsequently, a final 
requirement selects events for which the invariant mass 
of the two c-tagged jets is larger than 200 GeV. 

Simulated-event samples are used to aid the descrip¬ 
tion of the background and to model the SUSY signal. 
Top quark pair and single top production in the s and Wt 
channels are simulated with powheg-1.0 (r2092) [38], 
while the t channel single top production is simulated us¬ 
ing AcerMC 3.8 [39]. A top quark mass of 172.5 GeV is 
used. The parton shower, fragmentation, and hadroniza- 
tion are performed with pythia-6.426 [40]. Samples of 
IV-bjets, Z-bjets, and dibosons {WW, WZ, ZZ) with 
light and heavy flavor jets are generated with SHERPA 
1.4 [41], assuming massive 6/c quarks. Samples of Zti 
and Wti are generated with MadGraph-5.1.3.33 [42] 
interfaced to PYTHIA-6.426. The signal samples are gen¬ 
erated for a simplified SUSY model with only a sin¬ 
gle c state kinematically accessible, and with BR(c —>■ 
c -b Xi)=100%) using MadGraph-5.1.5.11 interfaced 
to PYTHIA-6.427 for the parton shower, fragmentation, 
and hadronization. Signal cross sections are calculated 
to next-to-leading order in the strong coupling constant, 
adding the resummation of soft gluon emission at next- 
to-leading-logarithm accuracy (NLO-bNLL) [43-45]. The 
uncertainty on each nominal cross section is defined by 
an envelope of predictions using different PDF sets and 
factorization and renormalization scales, as described in 
Ref. [46]. The Monte Garlo (MG) samples are processed 
through a detector simulation [47] based on GEANt4 [48]. 
The effects of pileup are included in the simulation. Ef¬ 
ficiency corrections derived from the data are applied to 
the simulation to correct for lepton efficiency as well as 
the tagging and mistagging rates. 

The main SM process contributing to the background 
after all signal region selections is .Z-bjets, followed by 
lU-bjets, top quark pair, and single top production. Most 
tt events contributing are tt —>■ bMvqq events, in which 
either a r lepton decays hadronically, or an e or /i is 
out of the geometric acceptance or not reconstructed or 
identified. Contributions from multijet, diboson, and as¬ 
sociated production of ti with W, Z are subdominant. 
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Noncollision backgrounds are found to be negligible. 

The estimation of the main background processes is 
carried out by defining a set of three data control regions 
(CR) that do not overlap with each other or with the sig¬ 
nal regions. The CRs are kinematically close to the SRs 
and each of them is enhanced in one or two of the back¬ 
grounds that is dominant in the SRs, while having low 
expected signal contamination (less than 1%). A statisti¬ 
cal model is set up in which the background expectation 
in the CRs and SRs depends on several parameters of in¬ 
terest: the normalizations of the dominant backgrounds, 
top {tt + single top), Z-l-jets and fT-|-jets, as well as on 
nuisance parameters including the effect of uncertainties 
on the jet energy scale (JES) and resolution, calorime¬ 
ter resolution for energy clusters not associated with any 
physics objects, energy scale and resolution of electrons 
and muons, c-tagging and mistagging rates, pileup, and 
luminosity. A profile likelihood fit of the background ex¬ 
pectation to the data is performed simultaneously in all 
CRs [49], and from it the background normalizations are 
extracted. The normalization factors, which are consis¬ 
tent with unity within uncertainties, are then applied to 
the MC expectation in the signal regions. 

The first control region is populated largely by tt and 
IT-l-jets. It contains events with exactly one isolated 
electron or muon with px above 50 GeV. The leading 
two jets, with px > 130 and 50 GeV respectively, must 
be c-tagged. To select events containing W —)■ ii', the 
transverse mass of the system is required to 

be between 40 and 100 GeV. The upper bound reduces 
possible signal contamination from SUSY models that 
produce leptons in cascade decays. Finally, it is required 
that A™'®® > 100 GeV and mcT > 150 GeV. The second 
control region is populated hy Z ^ events with 

two opposite-sign, same-flavor leptons, where the min¬ 
imum px requirement is 70 GeV for the leading lepton 
and 7(6) GeV for the subleading electron (muon). The 
transverse momenta of the leptons are added vectorially 
to the Axto mimic the Z —>■ i/D decay, and the mod¬ 
ulus of the resulting two-vector is required to be larger 
than 100 GeV. The leading two jets are required to be 
c-tagged and their px must each be above 50 GeV. The 
invariant mass mu of the two leptons is required to be 
between 75 and 105 GeV (Z-mass interval). A third con¬ 
trol region, populated almost exclusively by dileptonic tt 
events, contains events with two opposite-sign, different- 
flavor leptons, where the leading lepton has px > 25 GeV 
and the subleading lepton px is above 7(6) GeV for elec¬ 
trons (muons). It is required that > 50 GeV and 

mu > 50 GeV. The leading two jets are required to be 
c-tagged and have px > 50 GeV. In all GRs, events with 
additional lepton candidates beyond the required num¬ 
ber of signal leptons are vetoed using the same lepton 
requirements used to veto events in the SRs. 

The subdominant background contributions from di¬ 
bosons, Zti and Wtt are estimated by MC simulation. 


Finally, the residual multijet background is estimated 
using a data-driven technique based on the smearing of 
jets in a low-A™‘®® data sample with jet response func¬ 
tions [50]. 

The experimental and theoretical uncertainties affect¬ 
ing the main backgrounds are correlated between con¬ 
trol and signal regions, and the data observed in con¬ 
trol regions constrain the uncertainties on the expected 
yields in the signal regions. The residual uncertainty due 
to the theoretical modeling of the top-production back¬ 
ground is about 7%. It is evaluated using additional 
MC samples generated with AcerMC (where initial- 
and final-state radiation parameters are varied) an al¬ 
ternative fragmentation model (herwig), an alternative 
generator (mc@nlo), and by using diagram subtraction 
rather than diagram removal to account for the inter¬ 
ference between tt and single top kUt-channel produc¬ 
tion [51]. After the fit, the residual uncertainties on the 
lU-l-jets and Z-l-jets theoretical modeling account for less 
than 20% of the total uncertainty. The dominant contri¬ 
butions to the residual uncertainty on the total back¬ 
ground are from c-tagging (^20%), normalization un¬ 
certainties related to the numbers of events in the CRs 
(10%-20%), and JES (~10%). 

For the SUSY signal processes, theoretical uncertain¬ 
ties on the cross section due to the choice of renormaliza¬ 
tion and factorization scales and from PDFs are found 
to be between 14% and 16% for c masses between 100 
and 550 GeV. Prior to the fit, the detector-related uncer¬ 
tainties with largest impact on the signal event yields are 
those for c-tagging (typically 15%-30%) and JES (typi¬ 
cally 10%-30%). 

Table I reports the observed number of events and the 
SM predictions for each SR. The data are found to be be¬ 
low the SM background expectations, but consistent with 
them given the uncertainties. Figure 1 shows the mea¬ 
sured mcT and TOcc distributions in the toct > 150 GeV 
region compared to the SM predictions. Monte Garlo 
estimates are shown after the normalizations extracted 
from the profile likelihood fit are applied. For illustra¬ 
tive purposes, the distributions expected for the simpli¬ 
fied model with (c, Xi) masses of (400,200) GeV and 
(550, 50) GeV are also shown. 

Since no significant excesses are observed, the results 
are translated into 95% confidence-level (C.L.) upper lim¬ 
its on contributions from non-SM processes using the 
GLg prescription [52]. Figure 2 shows the observed 
and expected exclusion limits at 95% G.L. on the c-y? 
mass plane, assuming a single accessible c particle with 
BR(c —>- 0 - 1 - Xi) = 100%. The SR with the best ex¬ 
pected sensitivity at each point in the plot is adopted 
as the nominal result. In the region where the c-tagged 
analysis of the ATLAS i —>■ c + Xi search [18] provides a 
stronger expected limit, i.e. for mc — m^a < mw, that re¬ 
sult is used. The region excluded by the ATLAS monojet 
search described in Ref. [18] is shown separately as a grey 
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mcT (GeV) 

>150 

>200 

>250 

Top 

7.4 ±2.7 (7.1) 

3.9 ±1.6 (3.7) 

1.6 ±0.7 (1.5) 

Z-bjets 

14 ±3 (13) 

7.7 ±1.7 (7.0) 

4.3 ±1.2 (3.9) 

VF-l-jets 

7.2 ±4.5 (7.4) 

4.1 ±2.6 (4.2) 

1.9 ±1.2 (1.9) 

Multijets 

0.3 ±0.3 

0.2 ±0.2 

0.05 ±0.05 

Others 

0.5 ±0.3 

0.4 ±0.3 

0.4 ±0.3 

Total 

30 ± 6 

16 ± 3 

8.2 ± 1.9 

Data 

19 

11 
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TABLE I. Expected and observed number of events for an integrated luminosity of 20.3 fb“^ at = 8TeV. Top, Z+jets 
and fT+jets contributions are estimated using the fit described in the text. For comparison, the numbers obtained using MC 
simulations only are shown in parentheses. The row labeled “Others” reports subdominant electroweak backgrounds estimated 
from MC simulations. The total uncertainties are also reported. 



[GeV] 

FIG. 1. Distributions of mcT (top) and nice (bottom), and 
their corresponding SM predictions. Signal region selections 
(mcT > 150 GeV for the nice distribution) are applied, other 
than for the variable plotted. Arrows indicate the SR require¬ 
ments on mcT and nice- In the ratio plots, the grey bands 
correspond to the combined MC statistical and experimental 
systematic uncertainty. 


shaded area. Systematic uncertainties, other than in the 
c pair-production cross section, are treated as nuisance 
parameters and correlated when appropriate. For the 
SUSY scenario considered, the upper limit at 95% C.L. 
on the scalar-charm mass obtained in the most conserva¬ 
tive cross-section hypothesis is 540 GeV for m -o = 0 (in- 
creasing to 555 GeV for the central estimate of the signal 
cross section). Neutralino masses up to 200 GeV are sim¬ 
ilarly excluded for mg < 490 GeV. This significantly ex- 



FIG. 2. Exclusion limits at 95% C.L. in the c-y? mass plane. 
The observed (solid red line) and expected (dashed blue line) 
limits include all uncertainties except for the theoretical sig¬ 
nal cross-section uncertainty (PDF and scale). The band 
around the expected limits show ±1 (t uncertainties. The dot¬ 
ted lines around the observed limits represent the results ob¬ 
tained when moving the nominal signal cross section up or 
down by the ±1 ct theoretical uncertainty. 


tends the results of previous flavor-blind analyses [16, 17], 
which provide no exclusion for m^o > 160 GeV, nor for 
single light squarks with masses above 440 GeV. The sig¬ 
nal regions are used to set limits on the effective cross 
sections cr^jg of any non-SM processes, including the ef¬ 
fects of branching ratios, experimental acceptance, and 
efficiency, neglecting any possible contamination in the 
control regions. Values of larger than 0.44 fb, 0.36 fb, 
and 0.23 fb are excluded at 95% G.L. for mex greater 
than 150, 200, and 250 GeV respectively. 

In summary, this Letter reports results of a search for 
scalar-charm pair production in 8 TeV pp collisions at the 
LHG, based on 20.3 fb“^ of ATLAS data. The selected 
events have large A™'®® and two c-tagged jets. The results 
are in agreement with SM predictions for backgrounds 





























































5 


and translate into 95% C.L. upper limits on scalar-charm 
and neutralino masses in a simplified model with a single 
accessible c state for which the exclusive decay c —>■ c+Xi 
is assumed. For neutralino masses below 200 GeV, scalar- 
charm masses up to 490 GeV are excluded, significantly 
extending previous limits. 
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